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Radiative two-photon cascades from biexcitons in semiconductor quantum dots under resonant 
two-photon excitation are promising candidates for the generation of photon pairs. In this work, we 
propose a scheme to obtain two-photon emission that allows to operate under very intense driving 
fields. This approach relies on the Purcell enhancement of two-photon virtual transitions between 
states of the biexciton dressed by the laser. The richness provided by the biexcitonic level structure 
allows to reach a variety of regimes, from antibunched and bunched photon pairs with polarization 
orthogonal to the driving field, to polarization entangled two-photon emission. This evidences that 
the general paradigm of two-photon emission from a ladder of dressed states can find interesting, 
particular implementations in a variety of systems. 

PACS numbers: 42.50.Pq, 03.67.Bg, 78.67.Hc, 42.50.Dv 


INTRODUCTION 

The generation of non-classical states of light is a ma¬ 
jor goal in the implementation of photonic quantum tech¬ 
nologies mm- a case of particular interest is the gener¬ 
ation of photon pairs, since they present a wide range of 
applications in quantum information and quantum com¬ 
munications l 3 . Photon pairs are an important resource 
to generate heralded single photons [4] and are also used 
as a key element for quantum key distribution mm, 
quantum teleportation 0i or to implement entangle¬ 
ment swapping and quantum repeaters [MCE]. Numer¬ 
ous other examples, like quantum lithography [12], the 
absorption rate increase from organic molecules in two- 
photon microscopy nan], quantum walks of correlated 
photons [1J51 or the quantum computation of molecular 
properties m, illustrate the rich variety of applications 
that these non-classical states of light can find. 

Despite having such an impressive number of ap¬ 
plications, the number of ways in which these states 
can be generated is limited: most sources of photon 
pairs employed to date are based on parametric down- 
conversion m 0 CESHIEl- This mechanism can be imple¬ 
mented in several platforms, and thankfully for prospec¬ 
tive technologies, semiconductors are demonstrating ex¬ 
cellent performances HMD- These sources suffer how¬ 
ever from the major drawback that the number of photon 
pairs generated in each process shows Poissonian statis¬ 
tics, with a non-zero probability of having zero or more 
than one pair [22]. Promising candidates to overcome 
this difficulty are quantum emitters that naturally emit 
entangled photon pairs in a radiative cascade [23] . In 
the semiconductor case, the biexciton |B) in a quantum 
dot offers such an implementation, and emission of en¬ 
tangled photon pairs from the biexcitonic cascade has 
been demonstrated in recent years [241128] . As an alter¬ 
native to off-resonant excitation, it is possible to initialize 


the biexcitonic state by coherent two-photon excitation 
(TPE) [28ti32] . which increases the coherence and indis- 
tinguishability of the emitted photons as compared to 
non-resonant pumping. The generation efficiency and 
the indistinguishability of the photons can also be im¬ 
proved by bringing a cavity in resonance with the biex¬ 
citonic transition m to enhance the emission thanks to 
the Purcell effect [33] . A particularly interesting possi¬ 
bility is to place the cavity in resonance with half the en¬ 
ergy of the biexciton to enhance the rate of spontaneous 
two-photon emission, such that two photons are emitted 
simultaneously into the cavity mode 34 36 . The joint 
implementation of coherent excitation and Purcell en¬ 
hancement via cavity modes has already been discussed 
in the literature and shown to be promising [28, 35 . 

Under coherent excitation, the intensity of the pump¬ 
ing sets a limit to the repetition rate of two-photon gen¬ 
eration, since strong driving fields dress the excitons and 
spoil the biexcitonic structure m EH- On the other 
hand, a recent proposal [39] took advantage of such 
a dressing and demonstrated that a continuous source 
of TV-photon states—with photon pairs as the simplest 
realization—can be achieved in precisely this regime of 
strong admixing of the exciting laser with the emitter. 
This relies on a family of virtual two-photon transitions, 
so-called leapfrog processes [40] . in the strong driving of 
resonance fluorescence. Since virtual two-photon states 
are emitted away from the fluorescence peaks, they have 
a very small probability to occur on their own. These 
elusive photons are however precious [41, since they fea¬ 
ture giant quantum correlations and violate classical in¬ 
equalities [42] . Despite their scarcity, their existence has 
recently been demonstrated experimentally by measure¬ 
ments of frequency-resolved photon correlations [43]. It 
is therefore timely to capitalize on these photons to de¬ 
vise bright, continuous sources of photon pairs by har¬ 
vesting them with a cavity mode, with a Pur cell-effect 
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applied similarly to previous enhancements of quantum 
correlations [44j 05] from real photons emitted at the 
sidebands [4614481 . 

In this work, we bring together the three main ideas 
exposed above: i) TPE from the biexciton, ii) cavity 
Purcell-enhancement of virtual processes and iii) multi¬ 
photon emission from a dressed system. This realizes a 
versatile two-photon source operating in the continuous 
regime with a high repetition rate. In comparison with 
the case of a single dressed two-level system, the biexciton 
introduces an extra degree of freedom, the polarization, 
that provides a richer set of physical regimes. In partic¬ 
ular, we demonstrate the emission of degenerate photon 
pairs with polarization orthogonal to the laser—therefore 
suppressing the laser background and undesired excita¬ 
tion of the cavity — and different two-photon counting 
statistics, as well as emission of polarization-entangled 
photons. All these different regimes can be accessed opti¬ 
cally with the same sample just by changing the intensity 
and polarization of the excitation. This unprecedented 
versatility will push forward the generation and use of 
photon pairs in the laboratory. Even more importantly, 
it evidences that the fundamental concepts are suscepti¬ 
ble to be applied in different platforms, such as super¬ 
conducting circuits [49], and that new regimes of non- 
classical light emission are within reach with variations 
of the design. 

Our analysis starts with a general introduction of the 
model and follows with a detailed description of the fea¬ 
tures of the dressed biexciton alone, to finally move to 
the complete picture with the inclusion of a cavity that 
probes and enhances the single and two-photon transi¬ 
tions present in the dressed system. 


MODEL AND DRESSED STATE PICTURE 

The system under consideration is a semiconductor 
quantum dot with a biexcitonic structure, as depicted 
in Fig.[l|a). The Hamiltonian of this system is given by: 


H x = ( 1 ) 


where { 04 , 04 } are the annihilation operators of the ex- 
citons with spin {t,!}, is the excitonic energy (we 
consider degenerate excitons) and x i s th e biexcitonic 
binding energy. The biexciton frequency is, therefore, 
co>b = 2co>x — X* To separate the four-level system into 
two different polarization cascades we change to the lin¬ 
ear polarization basis: 


|H) = 7I (lt) + ll)) ’ |v> = 7l (lt) " lw 


with the annihilation operators 




<7 V = “ °l)- 


( 2 ) 


(3) 



FIG. 1: (Color online) (a) Biexcitonic level system in the lin¬ 
ear polarization basis (Horizontal-Vertical states). The two- 
photon laser excitation is represented with two curly blue ar¬ 
rows, at half the biexciton energy cjl = cjx —x/2. (b) Dressed 
state picture at strong laser pumping, where the vertical po¬ 
larization states (blue) transform into the new states |d=), |0)}. 
The three possible horizontally polarised transitions from 
these states to the |H), appear with curly arrows, (c) Spec¬ 
trum of emission in the two polarizations, horizontal (solid 
red) and vertical (dashed blue), and the three horizontally 
polarised transitions marked with vertical lines. Parameters: 
9 = 5 x 10 2 7 , x = 2 x 10 3 7 and g — 0. 


These operators describe transitions from the biexciton 
to an excitonic state or from an exciton to the ground 
state by emission of photons with the corresponding hor¬ 
izontal or vertical polarization (red and blue colors in 
Fig. 0 - We will neglect the small fine structure split¬ 
ting in frequency that is usually found between the two 
different excitonic states, since it has no impact in our 
scheme and only complicates the algebra. It can be triv¬ 
ially added if needed. 

We implement a continuous resonant excitation of this 
level structure that affects only one of the polarizations 
(chosen to be the vertical one without loss of generaliza¬ 
tion). This is accounted for by a coherent driving term 
in the Hamiltonian: 

H n = n (4e“ iWLt + °ve iULt ) (4) 

where Q represents the intensity of the exciting laser. In 
order to drive the biexciton state directly, the laser fre¬ 
quency is set at the two-photon resonance, cjl = cjb/ 2 . 
This results in a two-photon excitation (TPE) to the 
biexciton level [281432] . On the other hand, we gather 
and enhance the emission in the perpendicular polariza¬ 
tion (horizontal) through the coupling to a cavity mode 
with the same linear polarization. This way, we com¬ 
pletely separate in polarization the excitation and emis- 
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sion channels and do not need to worry about the elasti¬ 
cally scattered light from the laser. The coupling to the 
cavity mode is given by the Hamiltonian term: 

H c = (Jeo^a + ^(aVn + a cr^) (5) 

We write the total Hamiltonian in the rotating frame of 
the exciting laser: 

HmH 0 + H n + H c = 

Ax(0h0-h + d CT v) - X cr H <7 H 0 'v <7 V + fi(0y + O-y) 

+ Ac a 1 a + g{a ! cr-a + ao-jj) ( 6 ) 

where Ax = and Ac = uq— c^l- The dynamics of 

the whole system is described by a density matrix which 
follows the master equation: 

p=-i[H,p] + -C a p+'^ [A|X)<B| + £|G)<X|] p ( 7 ) 

z X=H,V 

where we use the definition of the Lindblad term: 

Cop = 2 OpO^ — O^Op — pO^O , ( 8 ) 

and the excitonic and cavity lifetimes are given by 7 and 
k respectively. We study the steady state of the system 
defined by p = 0 . 

Under TPE (Ax = x/2), the energy of the pho¬ 
tons from the laser matches half the biexciton energy, 
c.f. Fig. [lja). To understand the spectral features of the 
system before coupling it to the cavity (g = 0 ), we de¬ 
rive a dressed state picture for the biexciton [37 . The 
starting point is the set of bare states with n excita¬ 
tions, (|G)|n), |V)|n — 1), |H)|n — 1), \B)\n — 2 )}, where 
| n) describes the state of the driving field with n pho¬ 
tons. Since the laser is polarized in the vertical direc¬ 
tion, the state |H) is not dressed by it, while the rest of 
the excitonic states are. The new eigenstates are ob¬ 
tained by diagonalising the coupling Hamiltonian Hq 
(in the rotating frame of the laser) in the reduced ba¬ 
sis ||G)|n), |V)|n — 1 ), |B)|n — 2 )}, that is, the matrix: 

/ 0 n o\ 

h tpe = [n x /2 n . (9) 

\o n 0 / 

We do not include dissipation in this procedure since we 
consider it small as compared to U. This gives rise to the 
three new eigenvectors {|+), | 0 ), |—)} in each rung with 
the corresponding eigenenergies: 


A+ = ^ (-\/ x 2 + 32f2 2 + x) , 

(10a) 

A o = 0, 

(10b) 

A - = -1( 77+32^ - X ) , 

(10c) 


where the eigenvectors, dropping the photonic compo¬ 
nent from the notation, are given by |+) oc |G) + 
A+/fi|V> + |B), | 0 ) = (|B) - |V»/v/2 and |-) ex |G) + 
A_/U|V) + ]|B). Figure [ljb) depicts two successive rungs 
of excitation, including the state |H) which, as we said, 
remains bare. 

SINGLE PHOTON TRANSITION AND 
SPECTRUM 

The spectrum of emission of the system in each po¬ 
larization in the steady state, 5 x(cj), with X=H, V, is 
defined as Sx(w) = SR f™( a x(®) (J x(T))e lu,r dT. Both po¬ 
larizations are plotted in Fig. [ljc) for comparison. The 
number of peaks appearing and their positions can be 
explained in each polarization X by the allowed single 
photon transitions under the operator crx- In the case 
of H polarization, only transitions between |H) and the 
dressed states i = +,0, — are allowed: |<H|cth |^)| 2 7 ^ 0 
or | |o"h |H) | 2 7 ^ 0. The transition |H) |H) or be¬ 

tween dressed states |i) —>• | j) are forbidden in H po¬ 
larization, since |(H|<th|H)| = 0 and ||o"h| j)| = 0 for 
all = T, 0, —. The three possible transitions that 
can take place from the dressed states to |H), occur re¬ 
spectively at the following detunings from the laser (see 
Fig- 0b)): 

l+MIH): A I = i(y X 2 + 32fi 2 - X ) , ( 11 a) 

|0) —»• |H) : A„ =-x/2, (lib) 

|-)->|H): A III = -I(v'x 2 + 32n2+x) • (11c) 

The other three possible H-polarised transitions take 
place from |H) to the dressed states, at opposite detun¬ 
ings —Ai, —An and — Am. Remarkably, Su(uj) does not 
present any resonance at the laser energy. 

On the other hand, the spectrum in V polarization, 
Sy(a;), plotted with a dashed blue line in Fig. [ljc), con¬ 
tains seven peaks corresponding to the nine possible tran¬ 
sitions between dressed states, | i) —)► \j), with those three 
between the same dressed states, | i) -A |i), degenerate in 
energy at 

TWO-PHOTON TRANSITIONS AND 
SPECTRUM 

The next step in the characterization of the sys¬ 
tem is the calculation of the frequency-resolved sec¬ 
ond order correlation function or two-photon spectrum, 
#p 2 )(c<Ji, 0 J 2 ) [EffllSI], that conveys how likely is to detect 
two photons with frequencies uji, uj\ simultaneously. For 
that purpose we use a recently developed technique [SO] 
that makes the calculation of this quantity, previously 
very involved, computationally accessible, based on the 
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( a ) Frequency correlations U 2 ) 0 3 ) Cavity spectrum S(u) (a.u) 


(c) 



FIG. 2: (Color online) (a) Two-photon spectrum in H polarization for the TPE. In blue, sub-Poissonian statistics (antibunching), 
in red, super-Poissonian statistics (bunching) and in white, Poissonian statistics (uncorrelated). Parameters: x - 4 x 10 3 7, 
Q = 10 3 7 , T = IO 7 and g — 0. (b) Cavity spectrum of emission as a function of the cavity frequency u a in the strong coupling 
regime, g — 10 2 7 , k, = IO 7 . The plot on the right hand side shows the integrated signal, i.e., the cavity population n a . (c) 
Example of the two-photon transition |+) —X—X |+) in the H polarization. 


inclusion of the detectors in the system dynamics. The 
parameter T is the inverse response time of the detec¬ 
tor. It provides the frequency window in which photons 
are detected around cji, uj 2 . We fix it to an intermedi¬ 
ate value T = IO 7 , so that the detectors can resolve full 
spectral peaks (with width of the order of 7 ) without 
resulting in superimposed signals, 7 < T 12 . 

Figure [2|a) shows the H-polarized two-photon spec¬ 
trum from the light emitted by the dressed biexciton sys¬ 
tem. This map features seven antidiagonal red lines of 
super Poissonian correlations with g^ 1 (hyperbunch- 
ing) that correspond to a family of virtual two-photon 
processes that go from one state in a rung to another 
state two rungs below, jumping over any states from the 
rung in between (whence the denomination of leapfrog 
processes). It was recently demonstrated [42^ that this 
virtual character provides such strong quantum corre¬ 
lations that photon pairs can violate classical inequali¬ 
ties such as the Cauchy-Schwarz inequality. Whenever 
any two of the frequencies involved correspond to transi¬ 
tions between real states, these correlations change char¬ 
acter and the violation of Cauchy-Schwarz inequalities is 
spoiled. This can be seen in Fig. [2^a) as a piercing in 
the bunching lines whenever they intersect the vertical 
or horizontal ones, appearing at — cjl = d=Ai, ±An, 
±A m . 

Since the leapfrog lines originate from two-photon 
transitions, we can understand them in terms of the 
two-photon operator cthcth- Transitions starting or end¬ 


ing at |H) are not allowed, since |(H|<jh 0 'hK)| = 0 and 
|(^h^h|H)| = 0. All other nine two-photon transi¬ 
tions, | i) — X—X | j), occur, since 10 ’ |ct h cth |^) | 7 ^ 0 for all 
= T, 0, —, and give rise to seven lines which follow 
the general equation T (^2 — 2co>l = A 2P with: 

|i) \i) with z = +, 0 , — : A 2p = 0 , ( 12 a) 

|+) ^ |0> : A?x p = 1 (7x 2 + 32fi2 + , (12b) 

|0>->-H-): A X 2P = I(7F+32^-x) , ( 12 c) 

I+M-H-): *w = \Vx 2 + 3202. (! 2 d) 

The remaining three lines are described by inverting the 
order of the three last transitions and changing the sign 
of the corresponding A 2P . Figure J2^c) shows an example 
of a two-photon transition, |+) —X—X |+). 

In Fig. [3j we have another view of these leapfrog res¬ 
onances, selecting the diagonal of the two-photon spec¬ 
trum in Figure [2|a), that is, for c u = oq = u> 2 . The 
leapfrog processes appear as seven lines around 12 /x ~ 
10 “ 1 and spread as 12 is increased. The blue lines corre¬ 
spond to the single-photon resonances that are also ap¬ 
parent in the spectrum of emission, c.f. Fig. [3c). Re- 
ducing 12 below the dissipation levels (bottom part of 
the plot), the system experiences a transition into the 
spontaneous emission regime where there is no dressing 
of the levels and the spectral structures are much sim¬ 
pler: only two peaks for the spectrum of emission and a 
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single leapfrog peak at uo = cjl in the two-photon spec¬ 
trum. This regime was extensively investigated by one 
of the authors under incoherent excitation m- In the 
present work, where it appears as the low pumping limit, 
it will be used only for comparison with the high pumping 
regime. 


PURCELL ENHANCEMENT OF TWO-PHOTON 
TRANSITIONS BY A CAVITY MODE 

These virtual leapfrog transitions can be made real by 
coupling the system to a cavity (we switch on g ^ 0) 
in resonance with at least one of the two frequencies in¬ 
volved. If the coupling is strong as compared to the cavity 
dissipation, k, the two-photon emission can be Purcell- 
enchanced. We can observe this in the cavity spectrum 
of emission, given by S a (w) = / 0 °° (af (0)a(r))e luJT dr 
and plotted in Fig. |2|b). Because of the strong corre¬ 
lations between the two frequencies, the cavity Purcell- 
enhancement of one of the two photons of a bunching line 
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FIG. 3: (Color online) Two-photon spectrum with uj = uj\ = 
u 2 (the diagonal) of the biexciton system as a function of 
the driving field intensity Q. Top panel shows a cut along 
the dashed line in the bottom panel. Blue colors in the 
map represent sub-Poissonian statistics (antibunching), red, 
super-Poissonian statistics (bunching) and white, Poissonian 
statistics (uncorrelated). The physics changes from that of 
the biexciton spontaneous emission regime (bottom part), 
with a single leapfrog peak, to that of the dressed biexciton 
system (top part), with seven leapfrog peaks. Parameters: 
X = 4 x 10 3 7 , r — IO 7 and g — 0. 


triggers the emission of the second photon, even if this 
one is not in resonance with the cavity. This phenomenon 
leaves traces in the spectrum that help reconstruct the 
bunching lines when the spectrum is plotted as a func¬ 
tion of the cavity frequency. In this sense, the cavity is 
acting as one of the filters necessary to perform frequency 
correlations. 

As we have shown with coworkers in a recent work [39] , 
one can obtain useful two-photon emission by using this 
approach to Purcell-enhance two photons of the same fre¬ 
quency. This is evidenced by sharp peaks in the cavity 
population whenever it crosses one of the two-photon res¬ 
onances (Eqs. ( fl2| ) with uj a = = U 2 ), as can be seen 

in the right panel of Fig. |2|b). The single photon reso¬ 
nances appear as broad peaks and are detrimental for the 
two-photon emission. Therefore, the best candidates for 
pure two-photon emission are those leapfrogs far in en¬ 
ergy from other processes, that is, the sharp peaks with 
small overlap with the (one-photon) broad ones and that 
are further from other (two-photon) sharp ones. Logi¬ 
cally, it is also desirable that they are intense. The cen¬ 
tral peak, labeled I, at uj a = cjl is the best candidate 
since it is the most isolated one and is degenerate, with 
contributions from three different leapfrog processes. As 
we will discuss, this has consequences on the statistics of 
the emitted pairs. 

An accurate quantity to determine the quality of a 
two-photon resonance for two-photon emission is the pu¬ 
rity , 7T2 [39] . defined as the fraction of photons emitted 
in pairs from the total emission (including single pho¬ 
tons). Note that the purity being a probability, it is, 
unlike g^ 2 \ bounded: 0 < 7T2 < 1. Its definition is 
based on the fact that the photon counting distribu¬ 
tion of a perfect two-photon emitter shows a suppressed 
probability of counting an odd number of photons. The 
details of its definition and computation can be found 
in the Supplemental Material. In order to compute it, 
we simulate the actual emission of the system in the 
steady state via a quantum Monte-Carlo method [39] . 
The result is plotted in Fig. |4^d) for a cavity on reso¬ 
nance with each of the leapfrog peaks in the two-photon 
spectrum: I, II, III and IV, whose positions shift with 
U as plotted in panel (e). The corresponding cavity 
population n a = (a^a), second order correlation func¬ 
tion g^(O) = {af a 2 )/n 2 and two-photon second order 
correlation function [39] g^(0) = (aJ‘ 4 a 4 )/ (af 2 a 2 ) 2 ap¬ 
pear in (a), (b) and (c) respectively. The latter quantity 
takes the meaning of a standard second order correlation 
function for photon pairs when these pairs dominate the 
emission ( 7 ^ ~ 1). 

On the low driving regime, we can see that resonances 
I and II converge to the same point at uj = c^l and show 
very high purity: this is the usual regime of two-photon 
emission in the (undressed) biexciton, that has been stud¬ 
ied extensively before [52j[53]. Note, however, that this 
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high purity comes at the expense of the amount of signal 
(low n a ). As Fig. [4] shows, this signal can be enhanced 
by orders of magnitude if we increase the pumping in¬ 
tensity Q in order to bring the biexciton to the dressed 
regime. In this regime, all the resonances start being re¬ 
solved and the four of them present a sizable purity. In 
the case of resonances II, III and IV, the purity goes down 
whenever they cross a single-photon resonance (dashed, 
vertical lines in Fig[4|. At very high intensity, O > y, all 
of them reach almost 100% of pair emission. 

In this limit of high pumping, we observe a bunch¬ 
ing behaviour g^ 2 \0) > 1 for all the leapfrog resonances, 
which is an expected result for two-photon emission. The 



FIG. 4: (Color online) Steady state observables as a function 
of the pumping intensity Q for the parameters: y — 4 x 10 3 7 , 
g = 10 2 7 , k = IO 7 , A x = x / 2 an d A c = A 2P (blue), A^ 
(red), Ajn (yellow) and A 2 y (green). The gridlines mark 
the three points where leapfrog processes intersect with real 
transitions—dashed lines in (e)—, which spoils the purity. 


statistics of the photons, however, hides a non-classical 
behaviour if one regards the pairs as the basic entity of 
emission and consider the pair-pair coincidences as de- 
scribed by the g\ ; (0). In this case, we obtain antibunched 
photon pairs in resonances II, III and IV and bunched 
photon pairs in resonance I. This is not the only differ¬ 
ence between resonance I (blue) and the others. The 
position of this resonance is independent of the pumping 
intensity, and the emission at this frequency is order of 
magnitudes more intense than at the other resonances. 
These differences are explained by the fact that three 
different transitions contribute to the photon pair emis¬ 
sion at line I, and that the starting and ending state of 
the transition are always the same, as can be seen in 
Eq. (12a). Because of this, no reloading time to go back 
to the initial state is needed, which is the origin of the 
antibunching on all the other cases. All these features 
are a sample of the rich set of physical regimes that can 
be explored when the proposed method of multi-photon 
Purcell enhancement in dressed states systems is applied 
in non-trivial configurations. 


EMISSION OF ENTANGLED PHOTONS 


Many practical applications in quantum computing 
and quantum communication require emission of entan¬ 
gled photon pairs GSEffiU. So far, we have only consid¬ 
ered the case in which the emission was filtered by a sin¬ 
gle cavity with a fixed polarization. Therefore, changing 
the frequency of the cavity corresponds to moving along 
the diagonal of Fig. |2|a), and all photons extracted by 
the cavity will tend to be indistinguishable. However, 
the results for the spectrum of the cavity emission de¬ 
picted in Fig. [2|b) show that correlated photons of dif¬ 
ferent frequencies can be Purcell enhanced with a single 
cavity: therefore, two cavities in resonance with two dif¬ 
ferent, correlated frequencies—showing bunching in the 
map of Fig. <©( a)—will be expected to show strongly 
correlated emission. In our case, another interesting pos¬ 
sibility provided by the biexcitonic structure is to work 
with two degenerate polarizations of a single cavity, de¬ 
scribed by the bosonic annihilation operators an and 
ay. Two-photon emission then takes place in a reduced 
Hilbert space of polarizations {|HH), |HV), |VH), |VV)} 
of photons with the same energy. We will now show how 
the mechanism of two-photon emission described so far 
can also yield emission of entangled photons of the kind 
\ip) = (|HV) + |VH))/V2- The problem is theoretically 
described in the same way as before, but now with a 
different coupling Hamiltonian: 


Hq = CJc(«h«H + aytty) 

+ 5(4^ + «H °h) + S'fapv + ay cr^) ( 13 ) 
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The vertically polarized driving term used above leads 
to different probability of emission in horizontal or ver¬ 
tical polarization. Since we now want equal probability, 
we use a circularly-polarized driving laser: 

H n = n (ate- iWLt + a 0 e iWLt ) (14) 

with ctq = (<th + i(Jy)/y/2. We will not discuss in detail 
the possible single and two-photon transitions that arise 
in this model, since the physics and derivation are similar 
to the above material. However, it is now of interest to 
analyze the frequency-resolved, cross polarized second or- 

/Q\ 

der correlation function gf, (uq, UJ 2 ) of the dressed biex- 
citonic system alone (g = 0), which is a cross-correlation 
function between photons emitted at frequency uj\ and 
polarization H and photons emitted at frequencies UJ 2 



FIG. 5: (Color online) (a) Two photon spectrum for uj — 
uji — UJ 2 (diagonal) of the dressed biexciton for photons with 
opposed polarization (red) and same polarization (blue). For 
photons emitted at frequency A 2 y, the cross correlation be¬ 
tween polarizations is higher than the autocorrelations, (b) 
Autocorrelations and cross-correlations between two cavity 
modes with orthogonal polarizations coupled to the dressed 
biexciton as a function of their energy around A 2 y. (c) Con¬ 
currence of the emitted photon-pair state as a function of the 
total measurement time, (d) Density matrix of the emitted 
state for two different total measurement times denoted as n 
and t 2 . Parameters: Q = 8 x 10 3 7 , all the rest same as in 
Fig. [4] 


and polarization V. This correlation function can be com¬ 
pared with the frequency-resolved correlation functions 
for a fixed polarization that we have been considering so 

/'9'\ 

far, that we now term gf, ^ H (u; i,^)- Due to the circular 
polarized pumping, the system is symmetric under the 
exchange H H V, so ^rH H ( w i, ^ 2 ) = ^vv^ 1 ’^ 2 ) an( 4 

^r,Hv(^ij^ 2 ) = ^r,VH (kd ) ^2) • 

Figure [5|a) shows the comparison between these two 
quantities, with both photons having the same frequency 
(cji = CJ 2 = cj), equivalent to Fig. [ 3 } which was for 
the case of linearly polarized pumping. We observe that 
photons emitted at a frequency uj = cjl + A 2 y present 
strong cross-polarized correlations, clearly higher than 
the autocorrelation for each polarization. This is a non- 
classical result that violates the Cauchy-Schwarz inequal¬ 
ity [42], 54 . When the two degenerate modes described 
by an and ay are tuned to cjl + A 2 y and coupled to the 
system as described in Eq. (Jl3|, their cross correlation 

#hv = ( a H a v a H^v)/((^H^H)(^v^v)) and autocorrela¬ 
tion = ( a H a H a H&H)/( a H a H) 2 replicate the trend just 
described for the frequency-resolved second order corre¬ 
lations of the dressed biexciton alone. This can be seen 
in Figure [5jb), where the auto and cross-correlations of 
the two cavity modes are shown for a range of cavity 
frequencies around cjl + A 2 y. 

Quantum tomography [271 EDI S [56] allows us to re¬ 
construct the density matrix of the emitted photon pairs 
in the basis {|HH>, |HV), |VH), |VV)} from second or¬ 
der correlation functions corresponding to photon coinci¬ 
dence measurements. We define the (unnormalized) den¬ 
sity matrix: 


#AB,CD 



(A(0)apTOaD(T')ac(0))dT' 


(15) 


with A, B,C,D G {H,V}, where two-time correlation 
functions are calculated from the steady state of the sys¬ 
tem using the quantum regression theorem [57] . There¬ 
fore, r corresponds to the time of measurement that 
begins with the emission of the first photon, and for 
each r we define the normalized density matrix 9{r) = 
0(t) /Tt[0(t)]. This analysis reveals, for short measure¬ 
ment times, a highly pure density matrix, Tr[(9 2 ] ~ 0.92 
consisting of the entangled Bell state |^) = (|HV) + 
|VH))/V5 with fidelity F ~ 0.9, shown in Fig. |5[d). 
Beyond a certain time of measurement 1 /k, the density 
matrix loses purity due to the contributions from sub¬ 
sequent emissions. The degree of entanglement of this 
emitted bipartite state can be quantified by the concur¬ 
rence C ism, that in the case of pure states ranges from 0 
(separable states) to 1 (maximally entangled states) [58] . 
In our case, the concurrence takes a value C « 0.92 for 
short measurement times. However, one must bear in 
mind that the maximum concurrence for a mixed state 
is lower than one [59- and 0 is a mixed state with linear 
entropy Sl(0) = 4/3[l — Tr(# 2 )] ~ 0.11, which brings 

















this value of C closer to that of a maximally entangled 
mixed state. Not only is this result interesting by itself, 
but, being just a particular example, it also suggests that 
photons entangled not only in the polarization, but also 
in the energy degree of freedom—which in this case has 
been chosen equal for both photons for simplicity—could 
be obtained. 


CONCLUSIONS 

We have shown how Purcell enhancement of multi¬ 
photon resonances in the dressed ladder of a strongly 
driven biexciton can yield regimes of bright continuous 
two-photon emission. Thanks to the strong driving, the 
emission of photon pairs occurs at a much higher rate 
than it would in the approach that combines standard 
TPE (without dressing the system) and two-photon Pur¬ 
cell enhancement. The richness of the dressed biexcitonic 
structure allows to reach different two photon regimes 
like antibunched two-photon emission or entangled pho¬ 
ton pairs. These results suggest that the fundamental 
ideas behind this particular proposal are susceptible to 
be applied in a variety of platforms. 
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Enhanced two-photon emission from a dressed biexciton. 

Supplemental Material 

In this Supplemental Material we define the purity i r, a magnitude that quantifies the percentage of photons emitted 
in pairs from the total emission. This is therefore a bounded quantity i r G [0,1], and its definition is based on the 
observation that a perfect two-photon emitter will never emit an odd number of photons. We start the discussion by 
describing the probability distribution of the sum of two random variables in terms of generating functions. In our 
case, these two random variables correspond to the number of photons emitted by two different processes: one emits 
photons one by one, and the other, in pairs. The entirety of the light emitted is assumed to be the sum of these two 
processes, and we will show how this assumption does describes very accurately the photon counting distributions of 
the light emitted by the system. 


SUM OF TWO RANDOM PROCESSES 


For a given discrete random process X, we can define the generating function IIx = ($ X ), in such a way that the 
probability of obtaining X is given by P(X) = /ds x Ux\s=o- When we combine different random processes, the 

generating function is given by the product of the two original ones: Hx 1 +x 2 = Bx 1 Bx 2 - In the case of a coherent 
random process, the generating function is given by II a = e -A ( 1-s ). As we will see, this exponential form will be 
convenient. For the purpose of obtaining a general expression for the n-th derivative needed to compute P(X) from 
the generating function, we will use Faa di Bruno’s formula for the generalized chain rule expressed in term of Bell 
polynomials B n ^ k (x 1 , • • • ,x n - k +i) m - 


d n 
dx” 


f(g(x)) = '^2f {k \g(^)) B n,k ( g\x),g"(x ),• 
k= 1 


(n—/c+1) 


(SI) 


If /, as is the case for the coherent process, is just the exponential function, the n-th derivative of / appearing in (SI) 
can be taken out of the sum, which becomes just the sum of Bell polynomials known as the complete Bell polynomial 
B n (a i,--- ,a n ) = ^2k =1 Bn,k (&!,-•• ,a( n-/c+1 )). Therefore, for a generating function of the form II n = e 9 ^ s \ the 
corresponding probability for the random process is: 


1 d n e 0(°) 

P{n) = -r^-rn n | 5=0 = — —B n (ai, ■■■ , a n ) (S2) 

n! os n n\ 

with a n = g( n )(0), and F> 0 ({}) = 1- This reduces the problem of obtaining the photon counting distribution for a 
combination of random processes to express the generating function as a single exponential e 9 ^ and computing the 
n-th derivative of the exponent g^(0). 


THERMAL AND COHERENT DISTRIBUTIONS 


The approach to quantify the fraction of two-photon emission requires the selection of a proper ansatz for the 
photon counting distributions of single and two-photon processes, see Fig. [Si] For the case under discussion, we have 
obtained good results by considering a coherent one-photon distribution, a coherent two-photon distribution and a 
thermal two-photon distribution. The thermal component is essential to describe the photon counting distributions 
along the leapfrog line Aj P , giving incorrect results otherwise (see Fig. SI). The generating functions for a coherent 
one-photon process (Ai), a coherent two-photon process (A 2 ) and a thermal two-photon process (# 2 ) are given by: 


n Al = e“ Al(1 “ s) 

(S3) 

n A2 = e“ Al(1 “ s2) 

(S4) 

n 0 = l ~° 2 =e log (i- s 4 2 ) 

62 1 - s 2 e 2 

(S5) 

This give the total generating function II AljA2) 0 2 = II Al II A 2 IIe 2 = c: 9(s - ) . with 



1-02 


g(s ) = —Ai(l -s) - A 2 (l - s 2 ) + log 


1 - s 2 e 2 


(S6) 








iii 



FIG. SI: Photon counting distributions with the same purity according to some ansatz (in this case, coherent single and two 
photon processes) are sometimes badly fit. Upper row: Aj P . Lower row: Aj P . The case of higher Q for Aj P shows the ansatz 
is not correct in this case. 


whose n-th derivatives are: 


g( n )(Q) — 5 nt 1 Ai + S n 2 2A 2 + 


Q^ 2 e l[(n — 1)!] n even 


0 n odd 


This is all the information one needs to construct the photon counting probability given by Eq. (S2). 


(S7) 


THERMAL AND COHERENT COMPONENTS OF THE PURITY 


By fitting Monte Carlo photon counting curves such as the ones shown in Fig. [ST] one can obtain values for the 
parameters Ai, A 2 and # 2 for a given time window T. The mean values associated with each of the three processes 
are n\ x = Ai, n \ 2 = A 2 and Yiq 2 = 0 2 /(l — 0 2 ). The purity can then be defined as the fraction of the total mean value 
which is given by two-photon processes, that is: 


n\ 2 + n o 2 _ @2 / (1 ~ 02) + A 2 

■ no 2 Ai + A 2 + 62 /(1 — 62) 


n Xl + n \ 2 


This purity can be divided in a thermal plus a coherent part, tt = ttq + 7Ty, given by: 

02/(1 - 02 ) 


7T0 = 


7T0 = 


Ai + A 2 + 62/(1 — O 2 ) 
A 2 _ 

Ai + A 2 + 0 2 /(l — 0 2 ) 


(58) 

(59) 
(S10) 


The result of using the cothermal ansantz (S7) for the photon counting distribution and computing the purity as 


given by (S8) for the central leapfrog Aj in the biexciton configuration is summarized in Fig. S2 This provides the 


expected result of the purity being corresponding to that of a thermal process of bundle emission, which matches the 


/Q\ 

bunched values of g\ ; (0) in this line. Figure 
cothermal fit provides much better results. 


S3 


shows the associated errors in the fittings, demonstrating that the 
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FIG. S2: Purity for the central leapfrog as calculated by the standard method (7r) and by the cothermal ansantz ( 7 ^, 0 ), also 
with the separated thermal and coherent contributions. We see how the thermal contribution becomes important in this case 
and gives a higher purity than the one obtained with the standard, coherent ansatz. 
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FIG. S3: Error of the fitting procedure for the central leapfrog when using the standard coherent ansatz (blue) and the cothermal 
ansatz (red). This shows that the second method gives a much better fitting. 












